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ABSTRACT 

It is known that the class of luminous starburst galaxies tends to have higher R = 12 CO(J=l- 
0)/ 13 CO( J— 1-0) integrated line intensity ratios (R > 20) than normal spiral galaxies (R ~ 10). 
Since most previous studies investigated only R, it remains uncertain whether the luminous 
starburst galaxies are overabundant in 12 CO or underabundant in 13 CO. Here we propose a new 
observational test to examine this problem. Our new test is to compare far-infrared luminosities 
C/5 ; [^(FIR)] with those of 12 CO and 13 CO [L( 12 CO) and L( 13 CO), respectively]. It is shown that 

there is a very tight correlation between L( 12 CO) and L(FIR), as found in many previous 
studies. However, we find that the 13 CO luminosities of the high-i? galaxies are lower by a factor 
of three on the average than those expected from the correlation for the remaining galaxies with 
^ • ordinary R values. Therefore, we conclude that the observed high R values for the luminous 

, starburst galaxies are attributed to their lower 13 CO line intensities. 
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It is known that the class of luminous starburst galaxies with L(FIR)j> 10 n Lo tends to have higher 
R = 12 CO( J=1-0)/ 13 CO( J=l-0) integrated line intensity ratios than normal spiral galaxies (Aalto et al. 
1991, 1995, 1997; Casoli et al. 1991; Casoli, Dupraz, & Combes 1992a, 1992b; Hurt & Turner 1991; Turner 
& Hurt 1992; Garay, Mardones, & Mirabel 1993; Henkel & Mauersberger 1993; Henkel et al. 1998). It was 
often considered that the higher R values observed in these galaxies are attributed to the inflow of disk gas 
with high 12 C/ 13 C abundance ratios, possibly combined with a 12 C enhancement caused by nucleosynthesis 
in massive stars (e.g., Henkel et al. 1998). However, since most previous studies investigated only the ratios, 
it is uncertain whether or not these luminous starburst galaxies are overluminous in 12 CO or undcrluminous 
in 13 CO (Casoli et al. 1992a). In this Letter, comparing the luminosities of 12 CO(J=1-0), 13 CO(J=1-0) 
(hereafter 12 CO and 13 CO, respectively), and FIR emission for a sample of starburst and normal galaxies 
for which all the data are available in literature, we show that the observed higher R values are mainly 
attributed to a lower intensity in 13 CO with respect to 12 CO. We discuss possible physical mechanisms 
which can explain the 13 CO( J— 1-0) depression in luminous starburst galaxies. 



INTRODUCTION 



2. DATA AND RESULTS 

We have compiled 12 CO and 13 CO intensities from the literature (Wiklind & Henkel 1990; Aalto et 
al. 1991, 1995; Becker & Freudling 1991; Sage & Isbell 1991; Casoli et al. 1992b; Garay et al. 1993; 
Henkel, Whiteoak, & Mauersberger 1994). Our sample consists of 51 galaxies including ultraluminous 
infrared galaxies such as Arp 220. These data are used to estimate the CO luminosities; L(CO) is defined 
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as L(CO)= A x /(CO) K km s -1 pc 2 where A is the observed area in units of pc 2 and i(CO)= J Tpjj^ 1 ^ 
K km s _1 where Tj£ is the observed antenna temperature corrected for atmospheric extinction and -q is the 
main beam efficiency. 

The FIR data are compiled from the IRAS Faint Source Catalog (Moshir et al. 1992). The FIR 
luminosities are estimated using L(FIR)= 4ttD 2 1.26 x Kr n [2.58 x #(60) + 5(100)] (ergs s" 1 ) where 5(60) 
and 5(100) are the IRAS 60 /im and 100 fim fluxes in units of Jy and D is the distance of galaxies in units 
of cm (Helou, Soifer, & Rowan-Robinson 1985). Distances of nearby galaxies are taken from the Nearby 
Galaxies Catalog (Tully 1988); distances of other galaxies are estimated using a Hubble constant H = 75 
km s^ 1 Mpc -1 with I^gsr (radial velocity corrected to the Galactic Center: de Vaucouleurs et al. 1991). 
The compiled data are given in Table 1. Though our sample is not statistically complete, it is the largest 
sample compiled so far. 

As shown in Table 1, there are several galaxies which have R > 20. Hereafter we refer to these sources 
as high-i? galaxies (NGC 1614, NGC 3256, NGC 4194, NGC 6240, Arp 220, and Arp 299). ESO 541-IG 
23 [= Arp 236, R > 17 (lower limit)] and possibly IRAS 18293-3413 (i? ~ 26: Garay et al. 1993; R ~ 17: 
Aalto et al. 1995) are also high-i? galaxies. In total, there are eight high-i? galaxies in our sample. Except 
for the uncertain case of IRAS 18293—3413, all these high-i? galaxies are mergers. For the remaining 43 
galaxies, we obtain an average ratio, < R >= 11.3 ± 3.3. The values are consistent with previous estimates 
of 11 ± 3 (Aalto et al. 1991) and 9.3 ± 3.6 (Sage & Isbell 1991) for normal galaxies. 

First, we examine whether or not there is a beam size effect on R because the CO line observations did 
not cover the whole area of the individual galaxies. Ratios of beam size (Dbcam = half power beam width) 
to optical diameter at 25 mag arcsec -2 (-D25: de Vaucouleurs et al. 1991) are given in the third column of 
Table 1. In Figure 1, -Dbeam/-D25 is plotted as a function of R. We find no clear correlation, indicating that 
the difference in R cannot be attributed to the beam-size effect. 

In Figure 2, we compare L( 12 CO) with L( 13 CO). A tight correlation is found for galaxies 
with L( 12 CO)< 10 8 75 K km s _1 pc 2 . This almost linear correlation can be expressed as log 
L( 13 CO)=(0.96 ± 0.03) x log £( 12 CO) + (-0.68 ± 0.27). However, more luminous galaxies do not follow the 
same correlation, i.e., some of the luminous galaxies have lower 13 CO luminosities with respect to 12 CO. 

In Figure 3, we compare L(FIR) with both i( 12 CO) and L( 13 CO). It is shown that L( 12 CO) is well 
correlated with L(FIR); its correlation coefficient is 0.94. The average L( 12 CO)/L(FIR) for the high-i? 
galaxies, 0.013±0.007, is similar to that for the normal ones, 0.023±0.013. This correlation has been noted 
in many previous studies (e.g., Young & Scoville 1991 and references therein). In fact, Figure 3 shows that 
the correlation is established within a factor of three dispersion except for one galaxy (NGC 55) which is 
a nearby Magellanic- type irregular galaxy (Becker & Frcudling 1991). On the other hand, the correlation 
between L(FIR) and L( 13 CO) appears to be significantly poorer (the correlation coefficient is 0.86) than 
that between L(FIR) and L( 12 CO). The correlation is poorer because the high-i? luminous starburst 
galaxies have lower 13 CO luminosities than expected for galaxies with ordinary i? values; note that the 
average value of L( 13 CO)/L(FIR) is 0.00051 ± 0.00034 for the high-i? galaxies while that for the normal 
galaxies is 0.0023 ± 0.0017. Thus, we conclude that the higher i? in the luminous starburst mergers is 
mostly due not to enhanced 12 CO emission but to depressed 13 CO line intensities. An average i? value for 
the high-i? galaxies is ~ 30.7 ± 11.5. Since lower-limit data are used in this estimate, we obtain a 13 CO 
depression factor J> 3, given < i? >~ 11 for the normal galaxies. Any physical mechanism must explain the 
intrinsic weakness of 13 CO( J =1-0) emission by a factor of ^ 3 in the luminous starburst mergers. 
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3. DISCUSSION 

Here we investigate possible mechanisms that can explain the intrinsic weakness of 13 CO(J =1-0) 
emission by a factor of ;> 3 in the luminous starburst mergers. High R values may be caused by (see also 
Sakamoto et al. 1997): 1) low gas densities, 2) high gas kinetic temperatures, 3) large nuclear velocity 
gradients that keep CO column densities small at a given velocity, 4) small CO/H 2 abundance ratios, 5) 
more efficient excitation of 13 CO molecules toward upper levels (e.g., J > 2) with respect to 12 CO ones, or 
6) large 12 CO/ 13 CO abundance ratios. The first four mechanisms lead to lower opacity of the CO lines. 
The first mechanism also introduces non-LTE excitation, resulting in high R values. Since it is known that 
the luminous starburst mergers have a huge amount of dense gas (~ 1O 1O M ; Solomon et al. 1992; Scoville, 
Yun, & Bryant 1997), the first idea (lower gas densities) seems to be unlikely. Starburst galaxies often have 
both high gas kinetic temperatures and particularly large velocity gradients (Turner, Martin, & Ho 1990; 
Aalto et al. 1991, 1995; Devereux et al. 1994). However, if these mechanisms were responsible for the 
higher R values, we would find a 12 CO excess with respect to 13 CO. Therefore the second and third idea 
are also ruled out. The smaller CO/H 2 abundance ratio (i.e., lower metallicity) seems also to be unlikely 
for the luminous starburst mergers because the enhanced star-formation activity occurs in their nuclear 
regions where metal abundances are generally high (e.g., Storchi-Bcrgmann et al. 1996). It is known 
that 12 CO(J=3-2)/ 12 CO( J=l 0) ratios of starburst galaxies are often found to be higher than in normal 
galaxies (e.g., Devereux et al. 1994 and references therein). If 13 CO molecules in starburst mergers were 
excited more efficiently toward upper levels with respect to 12 CO, we would observe high-i? values in higher 
rotational transitions. However, there is no such evidence (Taniguchi, Ohyama, & Sanders 1998). Therefore, 
the remaining possibility is that the 12 CO/ 13 CO abundance ratios are systematically higher than those in 
normal galaxies. According to our result (Fig. 3), this should be achieved by a 13 CO underabundance with 
respect to 12 CO in the starburst mergers. 

We would like to thank Dave Sanders, Naomasa Nakai, Toshihiro Handa, Sumio Ishizuki, Naomi 
Hirano, and Satoki Matsushita for useful comments. We would also like to thank the referee, Christian 
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REFERENCES 

Aalto, S., Black, J. H., Johansson, L. E. B., & Booth, R. S. 1991, A&A, 249, 323 

Aalto, S., Booth, R. S., Black, J. H., & Johansson, L. E. B. 1995, A&A, 300, 369 

Aalto, S., Radford, S. J. E., Scoville, N. Z., & Sargent, A. I. 1997, ApJ, 475, L107 

Becker, R., & Freudling, W. 1991, A&A, 251, 454 

Casoli, F., Dupraz, C, & Combes, F. 1992a, A&A, 264, 49 

Casoli, F., Dupraz, C, & Combes, F. 1992b, A&A, 264, 55 

Casoli, F., Dupraz, C, Combes, F., & Kazes, I. 1991, A&A, 251, 1 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G. Jr., Buta, R. J., Paturel, G., & Fouque, P. 1991, 
Third Reference Catalogue of Bright Galaxies (Springer- Verlag) 



-4 - 



Devereux, N., Taniguchi, Y., Sanders, D. B., Nakai, N., & Young, J. S. 1994, AJ, 107, 2006 

Freeman, K. C, Karlsson, B., Lynga, G., Burrell, J. F., van Woerden, H., & Goss, R. 1977, A&A, 55, 445 

Garay, G., Mardoncs, D., & Mirabel, I. F. 1993, A&A, 277, 405 

Helou, G. I., Soifcr, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, L7 

Henkel, C, Chin, Y.-N., Mauersberger, R., & Whitcoak, J. B. 1998, A&A, 329, 443 

Henkel, C, & Mauersberger, R. 1993, A&A, 274, 730 

Henkel, C, Whiteoak, J. B., & Mauersberger, R. 1994, A&A, 284, 17 

Hurt, R. L., & Turner, J. L. 1991, ApJ, 377, 434 

Lonsdale, C. J., Helou, G., Good, J. C, & Rice, W. 1989, Cataloged Galaxies and Quasars Observed in the 
IRAS Survey, Version 2, JPL D-1932 (Pasadena: JPL) 

Moshir, M., et al. 1992, Explanatory Supplement to the IRAS Faint Source Survey, Version 2, JPL-D-10015 
8/92 (Pasadena: JPL) 

Rice, W., Lonsdale, C. J., Soifcr, B. T., Ncugebauer, G., Kopan, E. L., Lloyd, L. A., de Jong, T., & Habing, 
H. J. 1988, ApJS, 68, 91 

Sage, L. J., & Isbell, D. W. 1991, A&A, 247, 320 

Sakamoto, S., Handa, T., Sofuc, Y., Honma, M., & Sorai, K. 1997, ApJ, 475, 134 
Sanders, D. B., Egami, E., Lipari, S., Mirabel, I. F., & Soifcr, B. T. 1995, AJ, 110, 1993 
Scovillc, N. Z., Yun, M. S., & Bryant, P. M. 1997, ApJ, 484, 702 
Solomon, P. M., Downcs, D., & Radford, S. J. E. 1992, ApJ, 387, L55 

Storchi-Bergmann, T., Rodriguez-Ardila, A., Schmitt, H. R., Wilson, A. S., & Baldwin, J. A. 1996, ApJ, 
472, 83 

Strauss, M. A., Huchra, J. P., Davis, M., Yahil, A., Fisher, K. B., & Tonry, J. 1992, ApJS, 83, 29 

Taniguchi, Y., Ohyama, Y., & Sanders, D. B. 1998, in preparation 

Tully, R. B. 1988, Nearby Galaxies Catalog (Cambridge University Press) 

Turner, J. L., & Hurt, R. L. 1992, ApJ, 384, 72 

Turner, J. L., Martin, R. N., & Ho, P. T. P. 1990, ApJ, 351, 418 

Wiklind, T., & Henkel, C. 1990, A&A, 227, 394 

Young, J. S., & Scoville, N. Z. 1991, ARA&A, 29, 581 



This preprint was prepared with the AAS IAT^X macros v4.0. 



- 5 - 



Table 1. Properties of a sample of galaxies 
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Table 1 — Continued 
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Fig. 1. — Diagram of Dbeam/ D25 vs. R. Objects with R J> 20 (see the main text) are shown by filled circles 
and labeled. Note that IRAS18293 = IRAS 18293-3413 and ES0541 = ESO 541-IG 23. 
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Fig. 2. — Diagram of L( 12 CO) vs. L( 13 CO). The symbols have the same meaning as those in figure 1. A 
solid line corresponds to the correlation for R =< R > where < R > is the mean R using the data of galaxies 
with R < 20 (< R > = 11.3). The dotted lines corresponds to the correlation for R =< R > ±<r(i?) where 
a(R) = 3.3. 
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Fig. 3.— Diagrams of R (top), £( 12 C0) (middle), and L( 13 CO) (bottom) against L(FIR). The symbols 
have the same meaning as those in Figure 1. Alphabets in the top panel mean; a: NGC 4194, b: NGC 6240, 
c: NGC 3256, d: NGC 3256, e: NGC 3256, f: NGC 1614, g: NGC 3256, h: IRAS 18293-3413, i: Arp 220, 
j: Arp 299, k: Arp 299, 1: Arp 220, m: ESQ 541-IG 23, and n: IRAS 18293-3413. 



